Complex paths for regular-to-chaotic tunneling rates 
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In generic Hamiltonian systems tori of regular motion are dynamically separated from regions of 
chaotic motion in phase space. Quantum mechanically these phase-space regions are coupled by dy- 
namical tunneling. We introduce a semiclassical approach based on complex paths for the prediction 
of dynamical tunneling rates from regular tori to the chaotic region. This approach is demonstrated 
for the standard map giving excellent agreement with numerically determined tunneling rates. 
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Tunneling is a fundamental manifestation of quantum 
mechanics. Its basic features are established in standard 
textbooks [T], Q for particles confined by one-dimensional 
energy barriers: While the particle is classically trapped, 
it can escape quantum mechanically if the energy barrier 
is finite. This process typically exhibits an exponential 
decay exp(— 7t), which is characterized by the tunnel- 
ing rate 7. This rate describes the inverse life-time of 
the confined state and captures the relevant information 
of the tunneling process. Since time-independent one- 
dimensional systems have integrable dynamics, the tun- 
neling rates through the energy barrier can be computed 
from complex WKB-paths in the classically forbidden re- 
gion [iH3 
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Here, the imaginary part of the action S = J pdq, which 
increases with the width and the height of the barrier, 
is divided by Planck's constant h = 2TTti. In the semi- 
classical limit this ratio increases and tunneling vanishes 
exponentially. 

In contrast to one-dimensional potential wells, many 
systems of physical relevance have non-integrable Hamil- 
tonians such as driven atoms and molecules [H-Q, cold- 
atom systems Q , optical as well as microwave cavities 
and nano- wires [15|. These systems generically 
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have a mixed phase space in which classically disjoint re- 
gions of regular and chaotic motion coexist. Quantum 
mechanical transitions between such regions are called 
dynamical tunneling [l^, . In the above systems dy- 
namical tunneling is the key to understanding life-times 
and decay channels of long-lived states associated to a 
regular region. Furthermore, dynamical tunneling has 
important consequences for the structure of eigenfunc- 
tions [H, [1^ and spectral statistics [20M23| in mixed sys- 
tems. Due to this broad interest a lot of effort has been 
made to investigate dynamical tunneling experimentally 
[Mil and theoretically [M^. 

In this paper we derive a semiclassical prediction of 
dynamical tunneling rates 7 for the ubiquitous situa- 



tion of regular-to-chaotic tunneling. The focus is on the 
experimentally relevant regime, in which direct tunnel- 
ing to the chaotic region dominates. We generalize the 
WKB-formula, Eq. ([T]), to mixed systems by unifying 
the semiclassical time evolution method of complex paths 
[28j with the fictitious integrable system approach j36| . 
This demonstrates that direct regular-to-chaotic tunnel- 
ing rates are determined by complex paths, which con- 
nect regular tori to the boundary between the regular 
and the chaotic region, see Fig. [T] This approach is suc- 
cessfully applied to the standard map. 

In the following we explain our semiclassical complex 
path approach to dynamical tunneling rates for mixed 
systems. For simplicity we restrict the presentation to 




FIG. 1. (color online) Dominant complex paths (green ar- 
rows) for the direct regular-to-chaotic tunneling rate 71 of 
the standard map at «: = 2.9 and h = 1/50. The initial (in- 
ner red) torus and the final (outer blue) torus of the fictitious 
integrable system emerge from the real classical phase space. 
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time-periodic quantum systems, described by a time evo- 
lution operator U over one period of the driving. The cor- 
responding classical system is an area-preserving map U. 

We first introduce the fictitious integrable system ap- 
proach [3^. It is based on a fictitious integrable sys- 
tem H^cg{Q,p), which resembles the regular region of the 
mixed system as closely as possible and extends it be- 
yond its boundary, see Fig. [DJc). This allows to de- 
compose the Hilbert space into a regular and a chaotic 
part: As basis states of these parts we use the eigen- 
states of -ffrcgj which localize on tori with quantized ac- 
tion Im = 2^§p{<l)^1 = f^i'n^ + i)- The action /b of 
the boundary torus is used to divide these states into 
regular basis states \Im) with < Jb and chaotic ba- 
sis states with /,„ > /b, which will be called |/ch) in 
the following. Here, Ih is the (not necessarily quantized) 
action of the first torus of i/rog that is entirely located 
in the chaotic region of t/, see Fig. [2j In terms of the 
time evolution operator U of the mixed system, the rate 
7m of direct regular-to-chaotic tunneling from the mth 
quantizing torus to the chaotic region, is determined by 
summing the transition probabilities to all chaotic basis 
states [3™ 
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We now evaluate Eq. ^ semiclassically, which will 
lead to our main result Eq. (|6]). To this end we express 
the tunneling matrix elements (/ch|C^|^m) in the basis of 
the fictitious integrable system by generalizing the time 
evolution technique of complex paths [28j . Introducing 
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FIG. 2. (color online) (a) Phase space of the standard map 
at K = 2.9 with regular tori (lines) and a chaotic orbit (dots), 
(b) in the action-angle coordinates of -ffreg. (c) Phase space 
of -ffreg, (d) in action-angle coordinates. The boundary torus 
7b is marked by a thick line. 



the position representation gives 

(/ch|f/|/m) = / d<z' / dq {lM{q'\U\q){q\I,n). (3) 



Here, the first and the last factor are the initial and the fi- 
nal wave function of iJrog- These two factors are replaced 
by their semiclassical counter-parts, using the quantiza- 
tion of canonical transformations for the wave functions 
[soj . The propagator, {q'\U\q), is also expressed semi- 
classically [2^ . The arising integrals over oand q' 
are evaluated by the method of steepest descent [41|, i.e., 
a saddle-point approximation, which explicitly allows to 
take complex paths into account. This results in a semi- 
classical propagator for the tunneling matrix elements 
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It is constructed from classical paths i/, with action 
and a Maslov-phase shift 0^ . These paths v have to con- 
nect the initial torus /,„ to the final torus /ch • Since there 
are no such paths in real phase space, see Fig.[2l the above 
propagator is exclusively constructed from paths of the 
complexified phase space. Such a complex path v con- 
sists of three segments: (i) The first segment is the curve 
Cm,i/ on the analytic continuation of the initial torus Im 
of i/reg into the complexified phase space. This curve 
connects the real phase space to a specific point {q^,pi,) 
whose location is determined by the next segment, (ii) 
The second segment is the classical path of the complex- 
ified mixed system U. It has to connect {qu,Pu) on the 
initial complexified torus /,„ to a point {ql,p'^) on the fi- 
nal complexified torus /ch- This requirement determines 
the point {q^,p^) of the first segment, (iii) The final seg- 
ment is the curve Cch,^ which connects the point (q'^^p'^^) 
back to the real phase space, along the complexified final 
torus /ch of Hy-eg. The three segments are sketched in 
Fig. [1] (green curves). Their explicit determination for 
the standard map will be described below. 

The action S,^ of such complex paths is given by 

SAlch,Irn) = Jp{q)dq + S^{q,,ql) + jp{q)dq. (5) 



Here, the first and the last contribution arc action inte- 
grals carried out along the curves Cm.i/ and Cch,iy, respec- 
tively. The action SlJ{qij,ql) originates from the second 
segment, where the mixed system propagates between 
the initial and the final tori. From the possibly infinite 
number of paths, which contribute to the sum in Eq. Q, 
we select the dominant ones having the smallest positive 
imaginary action. To obtain the relative phase between 
the paths v, the curves Cm,i/ have to start at common ref- 
erence points and the curves Cch.i/ have to end at common 
reference points in real phase space. 
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In the following we evaluate the tunneling rates in 
Eq. ^ semiclassically. To this end, we replace the sum 
Sch integral Jj^ dich/fi, starting from the bound- 

ary action /b. Here, we exploit that transition matrix 
elements are semiclassically not restricted to quantizing 
actions. The matrix elements, Eq. are maximal at 
the boundary /b and decrease exponentially with increas- 
ing Ich- This allows to approximate the integral asymp- 
totically [4l|, leading to our final scmiclassical result for 
dynamical tunneling rates 
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Here, the only relevant paths are those which connect the 
regular action to the boundary action /b. Moreover, it 
generalizes the familiar WKB-prediction of Eq. (IT|) . The 
scmiclassical approximation leading to Eq. ^ cancels the 
interference terms between different paths. Hence, it is 
not necessary to determine the relative phases of the in- 
dividual paths for a scmiclassical prediction of dynamical 
tunneling rates. 

We illustrate our approach and demonstrate its suc- 
cessful predictions by applying it to the standard map 
(4^ . which is the paradigmatic example of a mixed sys- 
tem. It is obtained from the one-dimensional kicked 
Hamiltonian H{q,p) ^ T{p) + V{q)Y,^5{t - n) with 
T{p) = and V{q) = K/(27r)2 cos (27rg). The stro- 

boscopic time evolution is given by the map U 




FIG. 3. (color online) Dominant complex paths (green ar- 
rows) for the tunneling rate 71 of the standard map at k = 2.9 
and h = 1/50 as in Fig. [1] but in action-angle representation. 
The planes are the complexifications of the initial torus I\ and 
the final torus /b of -ffrog. White dots represent points on Ji 
which map onto /b. They are at the intersection of the gray 
lines (Im E' = 0) with the boundary between light red regions 
{ReE' < Eh) and dark blue regions (ReE' > E^)- The visu- 
alization is simpUfied by shifting the angle of the endpoints 
of step (ii) by —tuh along /b [as also done in Fig. [1]. 
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q + T'{p) p'=p-V'{q'). 



(7) 



The quantum time-evolution operator U is given by 
U ~ exp {—iV{q)/h) exp {—iT{p)/h), where h = 27rfi, is 
the effective Planck constant. We determine the tun- 
neling rate 7™ of the mth quantizing torus semiclas- 
sically, using Eq. One constructs a fictitious inte- 
grable system H^cgiq,p), Fig- c) , based on the frequen- 
cies of the main regular island [36|. For this H,-cg{q,p) 
there is a canonical transformation to action-angle coor- 
dinates (1,0), see Fig. [21d), giving H^cgil)- This allows 
for searching the complex paths v which connect the real 
initial torus Im and the real boundary torus /b of iJreg, 
according to steps (i)-(iii). For step (i) the torus Im is 
analytically continued by complexifying the angle 6. This 
torus forms a plane in the complexified phase space, see 
Fig. [21 parametrized by the real and the imaginary part 
of the angle 8. To obtain the canonical transformation to 
the corresponding (q,p)-coordinatcs, we numerically in- 
tegrate Hamiltons equations for H^cgiq,?) in imaginary 
time, starting from the real torus Im- The integration 
is done up to time t = ilm(0)/wm where cjm is the fre- 
quency of the real torus Im- For step (ii) the points 
{q,p) of the complexified initial torus Im are mapped 
by the complcxification of the standard map, Eq. ([7|), 
to points {q',p'). The points {q,p) are mapped onto the 
complexified torus /b, only if the energy E' = Hj-^giq' ,p') 



equals the energy of the boundary torus Ih- These nu- 
merically determined points are labeled {qy,Py) and are 
shown in Fig. [3] as white dots. For step (iii) we integrate 
Hamiltons equations for iJrog in negative imaginary time 
to get back to the real torus /b- Combining the three 
steps (i)-(iii) gives the complex paths v. Finally we eval- 
uate the action of these paths v according to Eq. Q in 
which the second contribution for the standard map is 
given by [H, [H] 

S^{q.,q'^)^^^^^^:^-V{q'^). (8) 

From the infinite number of paths v contributing to 
the sum in Eq. (|6]) we select the dominant ones, which 
have the smallest positive imaginary action, see arrows 
in Fig. [H There, the four paths v (together with their 
symmetry partners) contribute 20 %, 6 %, 33 %, and 
41 % (left to right) to the tunneling rate 71. Note, that 
we account for the parity of the standard map by just 
considering one symmetry partner and doubling its con- 
tribution in Eq. which leads to an additional factor 
of four in Eq. ([6]) . 

The semiclassically predicted tunneling rates 7,„, 
Eq. (|6]), for the standard map are shown as lines in Fig. m 
Numerically determined tunneling rates (dots) are ob- 
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FIG. 4. (color online) Tunneling rates 7™ of the standard map 
at K = 2.9 for m = 0,l,2,3vs. The semiclassical predic- 
tion, Eq. (|6|, (straight lines) is compared to numerical rates 
(dots) and for m = to the quantum evaluation of Eq. ([2]) 
(dashed line). The insets show the Husimi representation of 
the eigenstates with m = and m = 3at/i=l/60 and the 
corresponding quantizing tori (gray lines) in phase space with 
absorbing regions (gray). 



tained by opening the system [301 tangential to the reg- 
ular region (inset in Fig.|4]), which is consistent with the 
choice of /b (Fig. [5]) . The semiclassical prediction accu- 
rately describes the numerical tunneling rates, with de- 
viations by less than a factor of 2. This considerably im- 
proves the quality of previous predictions [36j , in which 
one finds quantization steps from the quantum mechan- 
ical evaluation of Eq. ([2]), sec Fig. SI These quantization 
steps do not occur in the semiclassical evaluation, as /b 
does not have to obey any quantization condition. The 
small deviations at small 7 can be attributed to non- 
linear resonances. One should be able take this effect 
into account, by taking our semiclassical prediction for 
direct regular-to-chaotic tunneling and combine [13] it 
with resonance-assisted tunneling. 

We now make a couple of remarks: (a) Our complex- 
path approach also works for other parameters k of the 
standard map, leading to predictions of similar qual- 
ity. Deviations due to nonlinear resonances may arise at 



larger 7 already, (b) Partial barriers which limit trans- 
port in the chaotic region can additionally suppress dy- 
namical tunneling [2^. This effect is not yet considered 
in our investigations since we open the system at the 
boundary of the regular region, (c) We obtain the tunnel- 
ing rates 7™ from semiclassical propagation over a single 
time step and find that only a small number of complex 
paths V contribute to Eq. ([5]) . Thus the determination of 
tunneling rates turns out to be considerably simpler than 
the semiclassical propagation of wave packets over long 
times, for which a careful selection among hundreds of 
domina nt p aths is needed [2^. (d) The natural bound- 
ary [4^ lU makes it impossible to continue regular tori of 
non-integrable systems deep into the complexified phase 
space [ij, [13, . We overcome this problem by com- 
plexifying the tori of the fictitious integrablc system iJrcg- 
(e) The choice of iJrcg is not unique and defining a pre- 
cise classical criterion for its quality is an open question. 
While the paths v may depend on this choice, we ex- 
pect to obtain a prediction of tunneling rates 7™ with 
similar agreement to numerical rates, (f) Determining 
tunneling rates just from WKB-paths along the initial 
torus Im of Jfreg; w hich was successful in specific situa- 
tions [10, m, [1J[1^ , gives for the standard map a rough 
estimate only. 

In summary, we have presented a complex-path ap- 
proach which allows for predicting regular-to-chaotic tun- 
neling rates in systems with a mixed phase space. We 
have successfully applied this method to predict tunnel- 
ing rates of the standard map. For the future it is de- 
sirable to develop a complete semiclassical description 
of rcgular-to-chaotic tunneling even in the presence of 
resonance-assisted tunneling. Moreover, the application 
of our approach to higher dimensional systems like bil- 
liards, optical microcavities, or atoms and molecules re- 
mains an open problem. Finally, we believe that our ap- 
proach is the basis to reveal the universal classical prop- 
erties which govern dynamical tunneling. 
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